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Abstract 

In this note we show that it is impossible to have data independent non-oscillatory three 
point finite difference scheme irrespective of its accuracy for a scalar hyperbolic initial value 
problem. 
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1 Introduction 

It is well elaborated in literature that evolution of discontinuities are inevitable in the solution of 
scalar hyperbolic conservation laws and artificial induced oscillations may appear in a numerical 
approximation in the vicinity of discontinuity EEllls]. One of the main emphasis in the numerical 
approximation of hyperbolic conservation laws is to construct schemes which are capable of yield¬ 
ing non-oscillatory approximation mum see also m- The numerical oscillations by high order 
schemes are well understood e.g., using modified equation analysis. In [T], Lax showed that induced 
oscillations analogous to the Gibbs phenomena must be present when the solution is approximated 
by a difference scheme that is more than first order accurate. Recently the cause of induced oscil¬ 
lations in the solution even by first order monotone finite difference scheme for hyperbolic problem 
is investigated in m 0 0 E]. This work starts with one motivational example to show that in¬ 
duced oscillations depend on the initial data. We define the notion of data dependent stability and 
uniformly non-oscillatory approximation. We finally show that it is impossible to have uniformly 
non-oscillatory approximation by any three point scheme irrespective of its accuracy. 


2 Uniformly non-oscillatory approximation 


We consider the following simple linear initial value problem, 


du{x, t) 
dt 


+ 0 , 


du{x, t) 
dx 


0, a 7 ^ 0, u{x, 0) = uo{x) 
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where u{x,t) is a scalar function of space x and time t variables, and uo{x) is a piece-wise smooth 
function with possible discontinuities, and characteristic function a{x, t) is a smooth function. We 
discretize the space and time variable in to a computational mesh using spatial and temporal step 
h and k respectively with grid points (xj = ih,tn = nk). Let A = ^ and uf represents and 
approximation of the value u{ih,nk). 

2.1 Lax-Wendroff (LxW) scheme 

Consider LxW scheme for Q. It is second order accurate therefore introduces oscillations in com¬ 
puted solution for discontinuous initial condition and is well understood by [T]. We now consider, 
Q corresponding to smooth initial condition 

a Uo{x) = sin(7rx), x G [—1 : 1] 

-1 

expi-^^ ifx G [—1 : 

0 else 

The numerical results by LxW are given in Figure It is clear from the results that even though 
both initial conditions are smooth, LxW does not introduces oscillations in solution corresponding to 
initial condition (a) while oscillations appears in the solution for initial data (b). Also as mentioned 
above even a hrst order monotone Lax-Friedrichs scheme exhibits local oscillations for discontinuous 
solution see [5]. Therefore it can be concluded that the induced oscillations depends on the initial 
data and can induce in to a numerical solution irrespective of its accuracy. This makes it reasonable 
to analyze the oscillatory behavior of any scheme with respect to the initial data. 





Figure 1: Solution using data CFL = 0.8, t = 6, N = 100. Both initial condition are with out 
discontinuity (a) No induce oscillation (b) induced oscillations 


2.2 A non-oscillatory condition for two point 

The idea is to use the following simple concept from geometry: Let each Pi, P 2 denotes two points 
in a hnite affine plane then there affine combination is dehned by 


P = Pi + 5{P2-Pi) (2) 

Geometrically, P represents a point on the line L which passes through Pi and P 2 , and if 0 < h < 1, 
P lies on the line segment joining point Pi and P 2 . Then following the method of characteristics and 
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affine combination (|^, an initial data independent non-oscillatory solution of Q can be ensured by 
an explicit two point consistent difference schemes scheme of the form 


11 ^+^ = 1 “ DAu^, if a{x,t) > 0, 

^ \ u'j+DA+u'j, ifa(x,t)<0, ^ 

provided 0 < D < 1. The coefficient D depends on a, X. A simple Von-Neumann stability analysis 
also shows that ^ is stable provided D < 1. We call scheme (|^ uniformly non-oscillatory (UNO). 


3 Non-oscillatory condition for three point scheme 


In this section we show that it is impossible to have a initial data independent uniformly non- 
oscillatory scheme with higher then two points. In [T], Lax pointed out the use of the following 
smoothness parameter to devise hybrid scheme which can avoid the Gibbs phenomena. It is dehned 
similar to the one in dziEaE] 


= 


qn,— 


A+uf 

A+< 

A_u7 


if a > 0 
if a < 0. 


(4) 


Using above affine combination we dehne 

Definition 3.1. Consider a consistent three point scheme written in the form 


n+i _ J '“j + 41(aA; 9i)A_u^ ifa{x, t) > 0 

I u'j + D{aX',6i)Aq.u'j ifa{x,t)<0, 


(5) 


where i G {j,j ± 1} and the parameter coefficient D depends on CFL number aX as well smoothness 
parameter 0 - The scheme 0 is uniformly non-oscillatory (UNO) if 


0 < DiOi) < 1, V G 


( 6 ) 


In case @ holds only for 6i G M\S where S C M, scheme 0 is called Data dependent stable (DDS). 
Note that three point scheme using a condition on 6 ensures that the updated value satisfies 
the following local maximum principle 

<_i < < <, ^fa > 0, 

Ui < ifa<0. 

3.1 Three point upwind linear schemes 

A generic numerical flux function for any linear three point upwind scheme for Q can be given by, 

_ J ati” -|- if a > 0 




+ f^u^+2 */ « < 0 - 

where for consistency the coefficients must satisfy a -|- /? = aQ 
consistent discretization requires i.e., F{u,u) = au 


( 7 ) 
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scheme 

Order 

a 

/3 

CFL condition 

DDS bound 

Two point upwind 

First 

a 

0 

0 < aA < 1 

UNO 

Three point upwind 

Second 

^ > 0 

§>0 

0 < aA < ^ 

^n,H- ■ 2—3aA o' 

Beam-Warming 

Second 

f(3- Aa) > 0 

f(l - Aa) > 0 

0 < aA < 1 


2-aX{3-a\) 3-a\ 
aA(l—aA) ’ 1—aA 


Beam-Warming 

Second 

f(3- Aa) > 0 

f(l - Aa) < 0 

1 < aA < 2 

et e 

2-aA(3-aA) 3-aA 
aA(aA—1) ’ aA—1 



Table 1: Data dependent stability bounds for three point upwind schemes for a > 0. Beam-Warming 
scheme changes its DDS interval with respect to CFL. 


For a > 0, the resulting three point conservative difference scheme using flux ([^ is 

■"+!=<-A(aA_nr-/3A_<_i) 


or 


u. 




A-u: 




A_u" 


A_u? 


By Definition 3.1 approximation ^ is DDS provided 

0 < A(a - < 1. 


( 8 ) 

(9) 


or 


which satisfies 


or 


-a 


a 


1 — Aa . 




1 - XCX Y} ^ n ^ 

— .f/5<o 


(10a) 

(10b) 


For a < 0, resulting scheme using flux in ([^ can be written as 

A+u 


= u] + X 


a — (3- 




A+< 


A+n" 


which by similar calculation is DDS provided 


cn - 1 ~ Aa , ^ 

1 — Aa _ a 


( 11 ) 


(12a) 

(12b) 


Note that DDS conditions (10) and (12) gives conditions on the initial data in terms of smoothness 


parameter 0 such that local oscillations does not introduced by a three point upwind schemes. In 
Table Data dependent stability region in terms of 9 is given for classical upwind schemes. 
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3.2 Three point centred linear schemes 


Consider the following generic consistent numerical flux function of three point centred linear scheme 
for 0 

(13) 




n 
i •> 


where again for consistency a + f3 = a. 

• a > 0 The resulting conservative approximation can be written as 

which using Q reduces to 


■-+i = nr-A[aA+nr + ^A-<] 


u. 


n+l 


= — \ (^aO'l’ + pj A_tf”, a > 0 


by Dehnition 3.1 (15) is non-oscillatory stable provided 

1-/3A 


- /? < aC’" < 


A 


it reduces to 


/3 1 - A/3 

- - < C’ < ^ , a > 0 

a Xa 

1 - A/3 


Aa 


<c 


< — —, a > 0 
a 


which on inversion gives condition for (14) to be DDS 

aA 


C’" 


qn,- 


E ^-OO,- 
E ( -oo, 


a 


13 

Xa 


U 


1 - A/3 


1 -/3A’ 

u 


oo , a > 0 


a 


a < 0. 


a < 0 scheme (14) can be written as, 

’*+1 =uf + x(a + /30"’+) A+<, o < 0 


which using similar calculation is DDS provided 


0? 


0? 


E -oo, - 


/3 


a 


U 


A/3 


E -oo. 


A/3 

1 — Aa 


1 — Aa 
U 


, oo ) , if /3 > 0 


f3 

-, oo 

a 


if /3 < 0 


(14) 

(15) 

(16) 

(17a) 

(17b) 

(18a) 

(18b) 

(19) 

(20a) 

(20b) 


The DDS conditions (18) and (20) classify the data type in terms of the smoothness parameter 9 
such that local oscillations does not get introduced by three point centred schemes. In Table 
data dependent stability reason is given for classical three point centred schemes which justihes the 
numerical oscillations by these schemes irrespective of their accuracy. In order to achieve uniformly 
non-oscillatory (UNO) approximation by an upwind or centred scheme one needs to choose coeffi¬ 
cients a, /3 in numerical flux 0 or ( |l3[ ) such that resulting scheme is DDS V E M. Bounds in 
(10) and (12) or (18) and (|20|) shows that scheme three point upwind or centred are not UNO and 


will always introduce local numerical oscillation except for the choice a = a, /3 = 0 in ([^ or the 
choice a = 0, a > 0 and /3 = 0, a<0in(13). This along with consistency requirement results in 
to classical first order upwind scheme. 
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scheme 

Order 

CFL 

a 

/3 

DDS bound 

Lax-Friedrichs 

First 

0 < aA < 1 

^{aX-l)<0 

^ (ttA + 1) )> 0 

0.’+G(oo, 1] U [(^“^,oo)] 

FTCS 

Second 

0 < aA < 1 

§>0 

§>0 

’+ G [ OO, 1] U 2-aA>Oo) 

Lax-Wendroff 

Second 

0 < aA < 1 

1(1 -aA) > 0 

f(l + aA) > 0 

e }+aA U 2+aA) 


Table 2: Data dependent stability bounds for three point centred schemes under CFL |a|A < 1 


3.2.1 Example: A hybrid second order UNO scheme 

The second order essentially non-oscillatory ENO scheme |in| for transport problem is UNO under 
the CFL condition 0 < Aa < 

Proof: Let a > 0, then the ENO second order reconstruction is 


= 


a 3a I I I 

+ yUi if \Ui-l - Ui\ < \Ui - Ui+l\ 


a a 
2^* + 


else 


The resulting ENO scheme for the transport problem Q can be written as 

■^+1 = ui- DA-Ui 


u. 


where 


D = 


aX 

aX 

aX 

aX 


1 A-ni_i 3 

"2 2 


1 / A+Uj 

2 1 A_Ui 

1 A+Uj I 

2 A-Ui 2 


if |A_aj| < |A+ttj| and |A_ttj_i| < |A_Uj| 
if |A_ttj| < |A+tij| and |A_tij_i| > |A_rti| 
^ |A_ai| > |A+ttj| and |A_ttj_i| < lA.-Uj] 
if |A_aj| > |A+ttj| and |A_ttj_i| > |A_rii| 


( 21 ) 


Note under the CFL condition 0<Aa<^,0<Zl<l, thus scheme is UNO. 


4 Conclusion 

An approach to find non-oscillatory condition on initial data in terms of smoothness parameter is 
demonstrated for three point schemes. It shows that it is impossible to have a data independent 
non-oscillatory three point scheme. In future it will be inter sting to extend this approach to analyze 
ENO schemes for their non-oscillatory condition and construction of hybrid high order UNO scheme. 
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